figS1_data = readtable('figS1.csv');

figure('position',[180,135,450,208+2/3]);

plot(figS1_data{:,1},figS1_data{:,2},...
    'o','MarkerSize',6,'linew',1,'Color',[0 0.447 0.741],...
    'MarkerEdgeColor',[0 0.447 0.741],'MarkerFaceColor',[0.0745,0.6235,1]);
set(gca,'fontsize',11,'position',[0.1300    0.2231    0.7750    0.7019],'xtick',-20:5:20);
hold all;
ylim([0,2.5]);
xlim(20*[-1,1]);
xlabel('Probe detuning \Delta [GHz]');
ylabel('Optical Depth {\itd}(\Delta)');
grid on;

plot(figS1_data{:,3},figS1_data{:,4},...
    '-','linew',1,'color',[0.8500    0.3250    0.0980]);
plot(repmat(figS1_data{:,3},1,8),figS1_data{:,5:12},...
    ':','linew',1,'color',[0.8510    0.7294    0.6784]);
